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In this paper, we develop the notion of the Hnear atomic quantum coupler. This device 
consists of two modes propagating into two waveguides, each of them includes a localized 
(— I , and/or a trapped atom. These waveguides are placed close enough to allow exchanging 

energy between them via evanescent waves. Each mode interacts with the atom in the same 



^ ' waveguide in the standard way, i.e. as the Jaynes-Cummings model (JCM), and with the 

^ ' atom-mode in the second waveguide via evanescent wave. We present the Hamiltonian for 

the system and deduce the exact form for the wavefunction. We investigate the atomic 
inversions and the second-order correlation function. In contrast to the conventional linear 
coupler, the atomic quantum coupler is able to generate nonclassical effects. The atomic 
inversions can exhibit long revival-collapse phenomenon as well as subsidiary revivals based 
on the competition among the switching mechanisms in the system. Finally, under certain 
conditions, the system can yield the results of the two-mode JCM. 

0^ ■ 

O ■ PACS numbers: 42.50.Dv,42.50.-p 

> 

I. INTRODUCTION 

Quantum directional coupler is a device composed of two (or more) waveguides placed close 



fl. 



The rate of flow of 



enough to allow exchanging energy between them via evanescent waves 
the exchanged energy can be controlled by the device design and the intensity of the input flux as 
well. The outgoing fields from the coupler can be examined in the standard ways to observe the 
nonclassical effects. Quite recently, this device has attracted much attention in the framework of the 
optics communication and quantum computing networks [2], which require data transmission and 
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a. 



ultra-high-speed data processing [3|]. Furthermore, the directional coupler has been experimentally 



implemented, e.g. in planar structures [4], dual optical fibres |5j| and certain organic polymers 



For more details related to the quantum properties of the fields in the directional couplers the 



reader can consult the review paper 



7|] and the references therein. 



The interaction between the radiation field and the matter (, i.e. atom), namely, Jaynes- 

n 

Cummings model (JCM) [q], is an important topic in the quantum optics and quantum information 
theories [91]. The simplest form of the JCM is the two-level atom interacting with the single-mode 
of the radiation field. The JCM is a rich source for the nonclassical effects, e.g. the revival-collapse 

n n 

phenomenon (RCP) [10], sub-Poissonian statistics and squeezing Furthermore, the JCM has 
3een experimentally implemented by various means, e.g. one-atom mazer [121, the NMR refocusing 

1 n n n 

13l |. a Rydberg atom in a superconducting cavity [1^, the trapped ion |15| and the micromaser jl6l |. 
Various extensions to the JCM have been reported including the two two-level atoms interacting 

nn 

with the radiation field(s) [171. Il8l]. 

The trapped atoms or molecules are promising systems for quantum information processing and 
communications 1^. They can serve as convenient and robust quantum memories for photons. 



20( 1 . The subject of coupling cold 



providing thereby an interface between static and flying qubits 
atoms to the radiation field sustained by an optical waveguide has already appeared in various 
contexts. For example, hollow optical glass fibers were used to guide atoms over long distances 



2ll |. especially, employing red detuned light field filling out the hollow core [22, l23|]. Substrate 



based atom waveguide can also be realized by using guided two-color evanescent light fields 



Moreover, the coupling of atomic dipoles to the evanescent field of tapered optical fibers has been 



demonstrated in 
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26l |. In this respect the optical nanofibers can manipulate and probe single- 



atom fluorescence. Moreover, it has been suggested that using a two-color evanescent light field 
around a subwavelength-diameter fiber traps and guides atoms. The optical fiber carries a red- 
detuned light and a blue-detuned light, with both modes far from resonance. When both input light 
fields are circularly polarized, a set of trapping minima of the total potential in the transverse plane 
appears as a ring around the fiber. This design allows confinement of atoms to a cylindrical shell 
around the fiber j2'jj]. Additionally, it has been shown that sub- wavelength diameter optical fibers 
can be used to detect, spectroscopically investigate, and mechanically manipulate extremely small 
samples of cold atoms. In particular, on resonance, as little as two atoms on average, coupled to the 
evanescent field surrounding the fiber, already absorbed 20 of the tot al p ower transmitted through 



the fiber. By optically trapping one or more atoms around such fibers 



28l | , it should become possible 



to deterministically couple the atoms to the guided fiber mode and to even mediate a coupling 
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between two simultaneously trapped atoms [29[. This leads to a number of applications, e.g., in the 
context of quantum information processing, high precision measurements, single-photon generation 
in optical fiber or EIT-based parametric four-wave mixing [s^ using a few atoms around optical 
nanofibers. Inspired by these facts we develop here the notion of the atomic quantum coupler 
(AQC), for which the interaction mechanisms inside the waveguides and between the waveguides 
depend on both the atomic and bosonic systems. These mechanisms are more complicated than 
those in the JCM, as we shall show shortly. For the AQC we show that the atomic inversions 
can exhibit long revival-collapse phenomenon as well as subsidiary-revival patterns based on the 
switching mechanisms in the system. Furthermore, under certain conditions, the system can give 
the results of the two-mode JCM. Also, the system is able to generate nonclassical effects. It is 
worth mentioning that the inclusion of one atom in one of the ports of the non-linear coupler has 



been considered in [3]|. Nevertheless, the solution of the equations of motion there is obtained by 
the rotation of axes, which does not give complete information on the system. 

We restrict the study in this paper to the development of the Hamiltonian model, its dynamical 
wavefunction and how does it work. These issues are discussed in section II. Additionally, in 
section III, we study two quantities, namely, the atomic inversions and the second-order correlation 
functions. 

II. MODEL FORMALISM AND ITS WAVEFUNCTION 

In this section we describe the linear directional atomic quantum coupler (AQC) and derive 



its wavefunction. A 
directional coupler 



so we discuss some basic differences between this device and the conventional 
Thus it is reasonable to shed some light on the linear directional coupler. 



1st 



which is described by the following Hamiltonian [7[: 

H ^ 

= ^Wjotaj -I- X{aia\ + a\a2), (1) 
i=i 

where oi (a|) and 62 {0^2) are the annihilation (creation) operators of the first and the second 
modes in the first and the second waveguides with the frequencies wi and 002] A is the coupling 
constant between the waveguides. Basically this device operates as a quantum switcher since it can 
switch the nonclassical effects as well as the intensities of the modes propagating inside one of the 



waveguides to the other 



321 ]. In other words, it can not generate nonclassical effects by itself. For 



some reason that will be clear shortly, we calculate the mean-photon numbers for the Hamiltonian 
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FIG. 1: Scheme of realization of the Hamiltonian (3). It is composed from two optical waveguides (yellow 
color). The circles in these waveguides denote the localized and/or trapped atoms. Mode 1 (2) pumped by, 
e.g., laser sources propagates along the first (second) waveguide and interacts with the first (second) atom 
via the coupling constant Ai (A2). The interaction between the first and the second waveguide occurs via 
the evanescent wave with the coupling constant A3. The outgoing fields from the coupler can be measured 
in the standard ways, e.g., using photon detectors. 

([1]) when the two modes are in the states |a,0). Thus we arrive at: 

{a\{T)ai{T)) = \a\^cos^{T), {al{T)a2{T)) = \a\^ sm^{T), (2) 

where T = Xt. These equations indicate strong switching mechanism in the linear coupler, where 
the intensity |ap in the first waveguide has been completely switched to the other one. Moreover, 
the mean-photon numbers can not exhibit the RCP. 

Now we are in a position to develop the AQC, which is the main object of the paper. The atomic 
coupler consists of two waveguides, each of which includes a localized and/or a trapped atom. The 
wavegudies are placed close enough to each other to allow interchanging energy between them. The 
two atoms (in the different waveguides) are located very adjacent to each other. In each waveguide 
one mode propagates along and interacts with the atom inside in a standard way as the JCM. The 
atom-mode in each waveguide interacts with the other one via the evanescent wave. The fields 
exited from the coupler can be examined as single or compound modes by means of homodyne 
detection to observe the squeezing of vacuum fluctuations, or by means of a set of photodetectors 
to measure photon antibunching and sub-Poissonian photon statistics in the standard ways. The 
scheme for the AQC is depicted in Fig. 1. From this figure and in the framework of the rotating 
wave approximation (RWA) the Hamiltonian describing the AQC can be expressed as: 
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J- = Ho + Hi, 

(3) 

where i^o and Hi are the free and the interaction parts of the Hamiltonian, cx^ and cij are the 
Pauh spin operators of the jth atom (j = 1,2); aj (aj) is the annihilation (creation) operator of 
the jth-mode with the frequency tOj and tOa is the atomic transition frequency (we consider that 
the frequencies of the two atoms are equal) and Ai (A2) is the atom-field coupling constant in the 
first (second) waveguide in the framework of the JCM. The derivation of the JCM Hamiltonian is 
well known, e.g. [33]. The interaction between the modes in the two waveguides occurs through the 
evanescent wave with the coupling constant A3. This term is the only one, which is conservative and 
can execute switching between the two waveguides. Thus it plays an essential role in the behavior 
of the AQC. We should stress that the switching mechanism occurs through the two JCMs (in the 
two waveguides) and can be obtained by applying the RWA in each individual waveguide. In other 
words, the quantity \^{ai(j^^^ 02^^^'^ + d|(T'^''d2(5"'^^) is nonconservative and hence it is cancelled 
out. Finally, the treatment of the switching mechanism in ([3]) is related to the notion of coupler, 
however, the existence of atoms in the waveguides has been taken into account. In ^ the treatment 
is considered only at the moment when the two fields interacting with atoms in the waveguides. 
Also when we treat the atoms (fields) classically the Hamiltonian ^ tends to that of the linear 
directional coupler (two-atom interaction). 

The interaction of two two-level atoms with the two modes has been considered in the optical 
cavity earlier [isl. [3^. 35 1. however, in the sense different from that presented above. For instance, 



as a sum of two separate Jaynes-Cummings Hamiltonians to investigate the entanglement 



3|] as 



well as the entang. 
localized qubits 



lement transfer from a bipartite continuous- variable (CV) system to a pair of 
35( 1 . Also, the quantum properties of the system of two two-level atoms interacting 
with the two nondegenerate cavity modes when the atoms and the field are initially in the atomic 
superposition states and the pair-coherent state has been investigated in [18i]. 

Next, we evaluate the wave function for the Hamiltonian We assume that the two modes 
and atoms are initially prepared in the coherent states |a,/3) and in the excited atomic states 
1 61,62), respectively. For resonance case 2uJa = uJi + uj2 one can easily prove that [Hq^Hi] = 0. 
Under these conditions, the dynamical wave function describing the system can be expressed as: 
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I ^(*)) = Yl Cn,m[Xi{t,n,m) \ 61,62, n,m) + X2{t,n,m) \ ei,g2,n,m+l) 

n,m=0 

+ X3{t,n,m) I gi,e2,n + l,m) + X4{t,n,m) \ gi,g2n + l,m + 1)] , ('^) 

where l^f) stands for atomic ground state. From the Schrodinger equation we obtain the foUowing 
system of differential equations: 



iXi{t, n, m) = \2\/ra + \X2{t, n, m) + Xi\/n + lX3{t, n, m), 

iX2{t, n, m) = A2\/w~+TXi(t, n, m) + A3 Y^(n + l)(m + l)X^{t, n, m) + \i\Jn+ 1X4(^5 n, m), 

(5) 

iX'^it, n, m) = Ai\/n + lXi{t, n, m) + \?j^J{n+ l){m + \)X2{t, n, m) + \2\Jm + 1X4(^5 n, m), 

1X4(1, n, m) = \i\Jn+ lX2{t, n, m) + \2\/ra + \Xj,{t, n, m), 

where the superscript "." means differentiation w.r.t. time. In the following, we give only the 
details related to the solution of the coefficient Xi{t,n,m), where the others can be similarly 
treated. Differentiating the first and last equations in ([5]) and re-substitute by the others we 
obtain: 

+ An,m)Xi{t,n,m) = -{i\'iC2D + ci)Xi{t,n,m), 
{0"^ + An,m)Xi{t,n,m) = -{i\3C2D + ci)Xi{t,n,m), (6) 
D = i,An,m = Af (n + 1) + Ai(m + 1), ci = 2X^X2^^ (n + l){m + 1), 02 = X^y^in + l){m + 1). 

From d^D one can easily obtain: 

{D^ + An,mfXi{t, n, m) = (iAgCzZ) + cifXi{t, n, m). (7) 

This equation can be easily solved. By means of the initial conditions stated above the exact forms 
of the coefficients Xj can be expressed as: 
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Xi(t,n,m) = ^exp{i^C2) cos{tQ-) — i-^?- sin{tQ-) +iexp(— i|c2) cos{t^l^) + i-^^ sin{tQ+) 



X2(t,n, m) 



2c2[A„,,„-4:^] 
^3 



33^{exp(iic2) [{4 - 2ci)Ai(m + 1) + (2^„,^ - ci) (A2C1 - ^(n + l)c2 



. sin(tn_ ) 

■ nz 



+ exp(-if C2) (ci + 2ci)Ai(m + 1) - (2^,^ + ci) A2C1 - ^(n + l)c2 



X3(t,n,m) 



-i\/n+l 



33^|exp(i|c2) (ci - 2ci)Af (n + 1) + {2An,m - ci) (AiCi 



sin(tr2+) 1 

—nr~ J' 



^(m + l)c2 



(8) 



, sin{tn_ ) 



X ' + exp(-i|c2) 



(ci + 2ci)Af (n + 1) - (2^„,„ + ci) ( Aici - A^(m + l)c2 



sin{tn-|- ) 



X4(t, n, m) = i exp(i|c2) — cos(tr2_) + sin(tJ7_) +iexp(— i|c2) cos(ir2+) + ^2^7- sin(tr2+) 



where 



1 



= -^A2(n + l)(m + l) + 4(AiV^H^±A2\A^H^)2. 



(9) 



It is obvious that the Rabi oscillation in the AQC is more complicated than that of the JCM. Prom 
the solution dS]) different limits can be checked. For instance, when (A2,A3) — > |0, 0) (A3 — > 0) 
the coefficients ([8]) reduce to those of the standard JCM (two decoupled JCM [34]). Moreover, 
when (Ai,A2) — > (0,0) the system reduces to a simple form, which is in a good correspondence 
with the conventional coupler ([1]). Nevertheless, the device, in this case, is a rich source for the 
nonclassical effects. This depends on the types of initial atomic states and can be explained as 
follows: (i) The atoms are initially prepared in |ei, 62). In this case the system reduces to the dark 
state, where i7mi|ei, 62) = 0. These states do not evolve in time. This property has been exploited 



in the quantum clock synchronization 



361 ]. (ii) The atoms are initially prepared in \e\^g'i)- The 



dynamical state of the system takes the form: 



*(r))= E Cn.rn cos[rv'(n + l)(m + 1)] I ei,52,n,m + l) 

n,m=0 '- 



(10) 



-i'sm\r^ (n + l)(m + 1)] | 5'i,e2,n + l,m) 
where T = Ast. The expression (|10|) reveals that the behavior of the radiation fields is typically 
that of the two-mode single-atom JCM [s^]. Finally, when the two atoms are initially in the Bell 
state [|ei,52) + \9\-,^'i)\l the wavefunction takes the form: 
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(b) 



0.15 
0.10 




FIG. 2: Evolution of the {cr'^P (t)) against the interaction time T = Xit with (a, /?) ~ (5, 5) for (A2, A3) = (1, 0) 
(a), (1,0.6) (b), (2,3) (c) and (1,1) (d). 

00 

l^(^)) = 75 E C„,mexpHrv^(n + l)(m + l)] [| ei, 92, n, m + 1)+ | 51, 62, n + 1, m)] . (n) 

n,m=0 

It is evident that the system exhibits atomic trapping, i.e. {a^z^\T)) = {a^\T)) = 0. Furthermore, 
the system is able to generate nonclassical effects, in particular, in the quantities, which depend 
on the off-diagonal elements of the density matrix such as squeezing (we have checked this fact). 

Now, we comment on the switching mechanism in the AQC. For the sake of comparison, we 
substitute /3 = in relations ([1])-(I8|) and calculate the mean-photon numbers as: 
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{a\{T)ai{T)) = \a\^ + ^ \Cnfl\^[\Xs{T,n,0)\^ + \X^{T,n, 



n=0 



(12) 



{al{T)a2{T)) = ^ |C7„,oR|X2(T, n, 0)|2 + |X4(T, n, 

n=0 



whre T = tAi. From these equations it is obvious that the intensity of the mode in the first 
waveguide cannot be switched to the other one. This is in a clear contrast with the hnear directional 
coupler (compare ([2]) and (jl2p l. This behavior is related to the nature of the atom- field interaction 
mechanism, which is close to the classic Lee model of quantum field theory. Moreover, this behavior 
is still valid even if the interaction between the modes and the atoms in the same waveguide is 
neglected, i.e. Ai = A2 = 0. In this case, expressions (jl2p exhibit the well-known RCP of the 



standard JCM [lO|. The final remark, AQC is able to switch the nonclassical effects from one 
waveguide to another based on the values of the interaction parameters. This is remarkable from 
(jl2p . where the mean-photon number in the second waveg uide (4(r)d2(r)) can exhibit RCP even 
though the second mode is initially in vacuum state. On the other hand, assume that the mode 
in the first waveguide is initially prepared in the even coherent state, which can exhibit squeezing, 
while the second mode is still in vacuum state. In this case, the density matrix of the second mode 

takes the form: 

00 

P2 = Y1 |C^2n,o|'{[|^i(T, 2n, 0)|2 + |X3(T, 2n, 0)\^]\0){0\ + [\X2{T, 2n, 0)\^ + iX^iT, 2n, 0)|2]|1)(1|}, 

n=0 

(13) 



where |C2nP is the photon-number distribution of the even coherent sate. From (jl3p . squeezing 
cannot be switched to the second mode. Nevertheless, if the second mode is prepared in the 
coherent state, it can exhibit squeezing. In this case, the source of the nonclassical effects could be 
the switching mechanism between the waveguides or the nature of the atom-field interaction. 

Now, we use above relations to investigate the atomic inversions and second-order correlation 
functions in the following section. For the sake of simplicity we consider a and /? to be real. 



III. ATOMIC INVERSIONS AND SECOND-ORDER CORRELATION FUNCTION 

Atomic inversion of the standard JCM is well known in quantum optics by exhibiting RCP. The 
RCP has a nonclassical origin and reflects the nature of the statistics of the radiation field. The 
evolution of the atomic inversion has been realized via, e.g., the one-atom mazer jl^ and using 
technique similar to that of the NMR refocusing ^]. In this section we investigate the behavior 
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FIG. 3: Evolution of the single-mode second-order correlation function as indicated against the interaction 
time r = Alt with (A2,A3) = (1,0) (a), (1,0.6) (b), (2,3) (c)-(d). 

of the AQC by studying the evolution of the atomic inversions and the second-order correlation 
functions. As the system includes two atoms we have two types of the atomic inversion, namely, 
single atomic inversion and total atomic inversion (cj^(T)) = ^[{a^\T)) + ((5"i^^(T))]. From (j4j) one 



can obtain the following expressions: 



{^r{T)) = E \Cn,M\X^{T,n,m)\^ + \X2iT,n,m)\^ - |X3(r,n,m)|2 - \Xi{T,n,m)\% 

n,m=0 

{af\T)) = g |C„,^|2[|Xi(r,n,m)|2 - |X2(r,n,m)|2 + \X^{T,n,m)\'' - \Xi{T,n,m)W (M) 

ra,m=0 
oo 

{^z{T))= E |C„,^P[|Xi(r,n,m)|2-|X4(r,n,m)|2]. 

n,m=0 
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As we mentioned in the preceding section the conventional directional coupler cannot exhibit RCP 
in the evolution of the mean-photon numbers. Nevertheless, the standard JCM can exhibit RCP 
provided that the photon-number distribution of the initial field has a smooth envelope. Similar 
conclusion has been reported to the two-atom single-mode JCM [l^. For the AQC we have found 
when a = P and Xj ^ the different types of the atomic inversions (I14p provide quite similar 
behaviors. It seems that the contributions of the coherence coefficients X2,X3 are comparable. 
Moreover, one can easily prove when A3 = and Ai = A2 the atomic inversions reduce to that of 
the standard JCM (see Fig. 2(a)). It is worth reminding that for the standard JCM the revival 
patterns occur in the atomic inversion over certain period of the interaction time afterward they 
interfere providing chaotic behavior. Additionally, the revival time is connected with the amplitude 
a through the relation T^- = 2TT\/n ~ 27r|a| [l^. We proceed, for Ai 7^ A2 they provide different 
forms of the revival patterns. Here we restrict the attention to the atomic inversion of the first 
atom (see Figs. 2(b)-(d) for the given values of the interaction parameters). We study three cases 
based on the relationship between the strength of the switching mechanisms in and between the 
waveguides, namely, A3 < Aj,A3 = Aj,A3 > Xj. Comparisons between Figs. 2(b)-(d) and Fig. 
2(a) are instructive. From Fig. 2(b) one can observe that the atomic inversion, after the zero 
and first revival patterns, exhibits long series of the subsidiary-revival patterns (see the inset in 
Fig. 2(b)). This behavior is completely different from that of the JCM. This indicates that the 
nonclassical effects generated by this device can sustain for an interaction time longer than that of 
the JCM. It is worth mentioning that the subsidiary-revival patterns have been observed for the 
JCM against the squeezed coherent state [ssl. This has been explained in relation to the photon- 
number distribution of the initial states. More illustratively, the photon-number distributions of 
the squeezed states exhibit many peaks structure, each of which gives its own revival patterns in 
the evolution of the atomic inversion. These patterns interfere with each other to produce these 
subsidiary-revival patterns. Nevertheless, for the system under consideration the occurrence of 
these patterns is related to the switching mechanism between the waveguides (compare Figs. 2(a) 
and (b)). This mechanism reflects itself in very complicated Rabi oscillations O-t as well as in the 
double summations in the atomic inversions formulae (I14p . Fig. 2(c) presents the case when the 
coupling constants are different. It is obvious that the RCP is still remarkable and the subsidiary 
revivals are smoothly washed out compared to those in Fig. 2(b). Generally, we have found when 
'^3 ^ Ai = A2 the atomic inversion exhibits long RCP (see Fig. 2(d)). Above information indicates 
that the switching mechanism between the waveguides plays an important role in the behavior 
of the AQC. Actually, we have found difficulties in giving mathematical treatment for the RCP 
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presented by the AQC since the Rabi oscillation is rather complicated. 

Now we draw the attention to the second-order correlation functions for the single-mode case, 
which is defined as: 

g) [t] = 1) J = 1,2, (15) 

(a](t)a,(t))2 

where g?^ (t) = for Poissonian statistics (standard case) , g?^ (t) < for sub-Poissonian statistics 
(nonclassical effects) and {t) > for super-Poissonian statistics (classical effects). The second- 
order correlation function can be measured by a set of two detectors, e.g. the standard Hanbury 
Brown-Twiss coincidence arrangement. For the system under consideration, this quantity is plotted 
in Figs. 3 for the given values of the interaction parameters. From these figures it is obvious that 
the AQC is able to generate long-lived sub-Piossonian effects, i.e. g\ (t) < 0. Furthermore, the 
basic features of the dynamics are still similar to those of the atomic inversion. Fig. 3(a) presents 
the well-known shape of the second-order function of the standard JCM. When the switching 
mechanism between the waveguides is involved the long RCP is dominant in the evolution of the 
g^^^ {t). Nevertheless, the shape of this phenomenon is quite different from that in the corresponding 
atomic inversion (compare Fig. 2(b) to Fig. 3(b)). For instance, the revival times in the two 
quantities are different. Also, the number of the subsidiary revivals in the atomic inversion is 

(2) (2) 

greater than that in the corresponding g^ (t). In contrast to the atomic inversions, g\ {t) and 
g^^ (t) can provide different behavior for the same values of the interaction parameters. This fact 
can be realized by comparing Fig. 3(c) to (d). 

In conclusion, in this paper we have developed, for the first time, the notion of the AQC. We have 
explained how does it work. Also we have derived the exact solution for the equations of motion. 
In contrast to the conventional coupler the AQC can generate nonclassical effects. Nevertheless, 
the switching mechanism in the former is more effective than that in the latter. Furthermore, the 
behavior of the AQC is sensitive to the types of the initial atomic states. We have shown that the 
system can give the results of the two-mode JCM under certain conditions. Additionally, we have 
discussed the evolution of the atomic inversions and second-order correlation functions. These two 
quantities can exhibit RCP, long RCP and long subsidiary-revival patterns based on the values of 
the coupling constants. Second-order correlation function can exhibit long-lived nonclassical effects. 
From the information given in the Introduction one can realize that the AQC is in the reach of the 
current technology. Also it may be of interest in the framework of quantum information. 
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